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This report summarizes some of the work performed under NASA GRANT No. 

NGR 36-004-065 to the University of Cincinnati during the period from August 1, 1^03 
jl9T^. The report was prepared by Professor L. H. Sobel, the Project Mana- 
ger » and by Mr. B. L. Agarwal, a Ph. D. candidate. Valuable contributions 
to this research effort were also made by Dr. J, Williams and Mi's. K. 

Hennessy of NA3A-LRC, and by Dr. T. Weller, NRC Fellow at NASA-LRC . 

A number of investigations pertaining to the buckling of cylindrical 
panels and complete cylindrical shells were performed. Some of the results 
obtained to date are presented in the following sections; 


A. OPTIMIZATION OF COMPOSITE STIFFENED CYLINDERS 

B. EFFECTS OF BOUNDARY CONDITIONS ON THE BUCKLING OF AXIALLY COMPRESSED 
CYLINDRICAL SHELLS 

C. LOAD INTRODUCTION TECHNIQUES FOR BORON INFILTRATED ALUMINUM PANELS 
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SECTION A 


OPTIMIZATION OF COMPOSITE 
STIFFENED CYLINDERS* 


■ * an expanded version of this section will be submitted for presenta- 
tion at the AIAA/ASME/SAE l6th structures, structural Dynamics, and Materials 
Conference to be held in Denver, Colo, on May 27-29, 1975. 



OPTIMIZATION OF COMPOSITE STIFFENED CYLINDERS 


!• INTRODUCTION : - An optimization study of composite stiffened cylinders 
is discussed in this section. The mathematical model for the buckling 
analysis has been coupled successfully with the optimization program AESOP 
(Ref, l). The buckling analysis is based on the use of so called "smeared 
theory" as used by Block, Card, and Mikulas (Ref. 2) for the buckling of 
stiffened orthotropic cylindrical shells. The equations used by Block, 
Card and Mikulas are modified to accomodate the laminated construction of 
the shell walls . 

2. DESIGN VARIABLES :- The loading, radius and length of the cylinder 
are assumed to be known parameters. An optimum solution then should give 
the value of cross-sectional dimensions and laminate orientations. These 
will be design variables. 

Figure 1 shows the optimized cylinder. It is assumed that stiffner 

spacing £ and ring spacing £ are unknown design variables. The skin 
s r 

of cylinder is allowed to have three different laminate orientations 
Otg» and which are assumed to be completely arbitrary. It Bhould be 

noted that skin is assumed orthotropic and each layer balanced. Hence, 
so far, there are a total of 8 design variables as shown in figure 1. 
Dimensions of the stringers and rings are discussed next. 

Figures 2(a) and (b) show a sketch of the rings and stringers, respec- 
tively. The rings and stiffners are assumed to behave as one- dimensional 
members. At this point the results previously obtained for the stiffened 
flat-plate (Ref. 3) are used to reduce the number of design variables. 
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Since the stringers and rings have similar characteristics, it will suffice 
to discuss only one of them. The stringers are composed of ±U5° laminates 
and 0° laminates as is shown in Figure 2(h). It is assumed that b g and 
h g are unknown design variables which decides the size of the stiffner. 
b^ s and b^ g are assumed to be *^ s and *8b s respectively. These 
values are based on previous results obtained for flat panels and this 
leads to 5 design variables for the stringers and, similarly, 5 for the 
rings. Hence, a total of 18 design variables are chosen as a starting 
point for the optimization work. 

* 

3. BUCKLING MODES :- Five different types of buckling modes are considered. 
These modes are as follows: 

(a) Gross buckling 

(b) Panel buckling (buckling between rings) 

(c) Skin buckling (buckling of the skin between contiguous rings and 
stringers ) 

(d) Local buckling of stringers 

(e) Local buckling of rings 

The buckling loads are determined from the analysis given in the next section. 
k. THEORITICAL ANALYSIS 

The notations and sign convention used in Ref. 2 is employed herein. 
h.l CONSTITUTIVE EQUATIONS 

For a laminated shell the stress strain relations for p th layer are 


given by 
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then, for a symmetric laminate layup, it can be shown for a Donnell-type 
analysis that (see ref. it) 
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In these equations v* v, w and N , N , N , M , M , M are* 

x y xy x’ y xy 

incremental displacements and stress resultants that take place 
ling. There positive sign conventions are shown in Fig. 4. A 
partial derivative with respect to the indicated variable, and 


respectively, 
during buck- 
comma denotes 


A., = 


E (QjX (a - d .) 
- ij p p p-i 


p=i 


‘ £ 'V, 


a 3 ,) 

P-1 


i,J = 1,2 and 6, d^ is the distance of the center of p layer from 
the reference axis, and P is the total number of layers. 

In general d i 6 and terms are not zero, but in the present work 

they will be assumed to be zero. 


4.2 BUCKLING ANALYSIS 


* Gross, Panel and Skin Buckling 

With the above more general constitutive relationships, the buckling 
equation of Ref. 2 have been modified to account for the laminated wall 
construction. This yields: 



( 1 ) 
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In this equation, N and N are prebuckling stress resultants (from 

x y 

now on, N will be taken equal to zero), m is the number of axial half 

y 

waves , n is the number of circumferential full waves , and 



In these expressions, EA is the extensional stiffness of the stiffners, 

GJ is its torsional stiffness, and El is the bending stiffness of the 
stiffner about the skin reference surface. Subscript s and r, respectively 
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represent stiffner and ring. For the gross ‘buckling mode, the above equa- 
tion is used directly, but for panel and skin buckling modes it is modified 
slightly. For panel buckling, the length of the cylinder is assumed to be 
equal to ring spacing and all the ring stiffness properties are set equal 
to zero. And for skin buckling, all the terms due to stiffner stiffness 
and ring stiffnesses are set equal to zero and a buckling load corresponding 
to - 


L = £r 
n — Integer 



n = 1, 2, 


gives the skin buckling load. 

• Calculation of stringer and ring stiffness properties : 

Equation 1 requires knowledge of the stiffness properties of the rings 
and stringers. These will now be determined. Since the stringer and ring 
are similar geometrically, it will suffice to discuss the stringers only. 
Figure 5 shows three members of the stiffner, each having a width b^. 

The width of each element is given by 
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If [A], is the extensional stiffness matrix for i member then 

l 

xt. 

Young's Modulus E . . for i member in x direction is given by (see ref. b) 
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where h^ is the thickness of the i"^ member given by 


h = t_ + t 
1 Is 2s 


h 2 = t 2s 


h 3 = *26 + *38 


Let EA denote the extensional stiffness of each stiffner. Then the 
is 

total extensional stiffness of the stiffner. EA is 

* s 
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The distance of the stiffner neutral axis from the skin reference 
axis , z g , is given by 

- 2^*3+ (h S + V* h S 

Z s EA 

s 


The bending stiffness (El ) about the skin reference axis is given 

o s 


by 
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The torsional stiffness (GJ)^ of the stiffner is computed as follows 


{GJ) s 4 ( • 9b . h ) 2 L 
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The contribution of the inplane shear stiffness due to stiffner 
(Agg) s is given by 
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With the help of the above stiffness properties the gross, panel and 
skin buckling loads can be computed. 

In order to determine the local buckling of stringer and ring it will 
again be sufficient to discuss only one of them. 

• Local Buckling 

For the local buckling modes of the stiffners , the buckling of members 
2, 3 and the skin between the webs is considered. All these members are 
assumed to be orthotropic plate members simply supported on all four edges. 
Hence the buckling load for i member is given by (Ref. 5*) 

KT^ +d 12 + 2D 66], 

k.3 CONSTRAINT CONDITIONS 

• Buckling Constraints 

For an optimum design to be a valid design, the applied load carried 
by each member cannot exceed the buckling load of the corresponding member. 
These buckling loads are now computed from a membrane prebuckling deforma- 
tions , the relation between loads and strains can be written as 
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Hence 


{e } = [A]" 1 {N} 

P P 

The prebuckling strains can thus he computed from the total prebuckling 

stress resultants ' N and N . Next the loads carried by the skin and the 

xp yp 

individual members of the stiffners will be computed. 

• * Load carried by the skin 
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* » Load carried by stringer members 
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• * Load carried by ring members 
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• Material Failure Constraints 

The skin and stiffner laminates must be checked for possible material 
failure. For the case of laminated composite stiffened cylindrical shells 
it will be necessary to check the strain in each laminate for failure. This 
is the most conservative failure criterian and is used for the present 



11 


problem in view of the lack of any other presently satisfactory failure 
criteria. 

If the laminate fibers are oriented at an angle 0 from the axial 
direction, then the strains in that laminate are given by {see. Ref. k) 
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where is strain along the fiber, e strain perpendicular to the fiber 

and y ttLe siiear strain. The strains given by above equation are 

constrained to satisfy the yield strains of the material in each laminate. 

5. NUMERICAL RESULTS 

5.1 COMPUTER PROGRAM 

A priliminary listing of the computer program developed for the 
optimization studies is given in appendix A. Detailed documentation of the 
use of this program and its capabilities will be given at the completion of 
this continuing effort. 

5.2 CHECK CASES 

In the use of smeared theory it is required to calculate the extensional, 
torsional and bending stiffness properties of the stiffners. In order to 
assess the effects of modifying the equations of Ref. 2 and of the assumptions 
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made in calculation of the stiffness properties, the following three check 
cases were used: 

1. Unstiffened cylinder using BUCLAP 2 (Ref. 6) 

2. Stiffened cylinder using BUCLAP 2 

3* Stiffened cylinder using BUCLASP 2 (Ref. j) 

These cases are briefly discussed next. 

The first two check cases are used primarily to check the effect of 
modifying the equations of Ref. 2. BUCLAP 2 can be used for computing 
buckling loads of a sitffened cylindrical shell by adjusting the stiffness 
matrices to account for the effect of eccentricity and stiffners. These 
modified stiffness matrices, which were computed in a related panel buck- 
ling study, are given in appendix B. Buckling loads for both stiffened 
and unstiffened cylinder were found to be in good agreement. 

The third case is used to check the effect of the use of smeared 
stiffners instead of discrete sitffners, and also to check the assumptions 
employed in the calculation of the stiffness properties of the stiffners. 

It was disappointing to find that discrete theory using BUCLASP 2 
gave a buckling load 30$ lower than that predicted by smeared theory. The 
reason for this difference is probably in the computation of the torsional 
stiffness of the stiffner assumed for smeared theory. This contention is 
supported by the fact that, when the cylinder was forced to buckle in 
axi symmetric mode the buckling loads given by smeared and discrete theory 
were almost the same. An investigation is presently under way to resolve 
this problem. The optimization results for stiffened cylinders will then 
follow. However optimization studies for unstiffened composite cylinder have 
been successful. The results of these studies are given next. 
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5-3 RESULTS FOR UNSTIFFEEED CYLINDERS 

Some preliminary results for unstiffened composite cylinder are pre- 
sented in Figs. 6 end 7* 

Figure 6 shows the weight strength plot for unstiffened cylindrical 
shells under uniform axial compression. There is clearly a weight saving 
of about using Graphite/Epoxy over aluminum. Furthermore, a similar 
or even better weight savings can be expected for stiffened cylinder, be- 
cause Graphite/Epoxy stiffners can carry the load more effectively than 
aluminum stiffners. 

Figure 7 shows results for unstiffened cylinders with the material 
properties used by Dow and Rosen (Ref. 8). They showed an "isotropic" 
arrangement of fibers was most optimum. But the present results show that 
a more optimum fiber orientation can be obtained using a general fiber orien- 
tations. Dow and Rosen also showed that isotropic configuration was better 
than ± 15 ° configuration, but this was not found to be so in the present 
computations. Both ±15° and isotropic configuration gave the same value 
of weight parameter as indicated in Fig. 7» 
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APPENDIX A 


LISTING OF COMPUTER PROGRAM 




S UB UT I N E CVL QPT : 

REAl. ALPHA! 201 - — 

REAL FUNCTNU00.) 

£ yUJV ALE M C EJ. A OAT At 74 2? » ALPHAj 

E QU I VAL ENCE ( AO AT AJ_2736 1 , JJ Ji 

E QU I V AL ENCE t A D A TA1 2948 I #M AXJ JJ A 

' EQU I V AL E N C E (A D A T A ( 2 5 5 2 ) f „FU NCTNJ 

C OM MON /AESOP D / _ A DA T A l 5 0 0 0J_ 

c 0 M MO N / S U B / FT OL {.25 i r PS I rt T J. 2 5JL»_ S IBA R ( 2 5_)_ iTOL FACJ. 2_5 ) » T_T 0 L12 5>_j _ 

~~T WT DQ W Hi 2 5 J_, UTU PI25) »NPSi(25) ; 

COMMON /AREA/ TS t A$S,A$R 

' 0 1 M E N S 10 N. QR l 4 Lt PI 5 j. 4 L* AJLiij. J£ A5_L21j^ARi. 3J j£XSJL3JLt-RXRLL^ij-6LU.CKS.i.9.Lt^l — 

____i C R (3 ) r MCR ( 3 i t NCOUNT UJ» E P SL±A52j_EPSL2l5A 

~ COMMON /DAT A/ BUCKL « EA $ t EAR t B2$ 1 82R 

COMMON /Q/ QL (4,5) 

COMMON /MAT/ EX » E Y ♦ ENUXV , GXY »RQ ______ — 

COM MON/ SAVE U E AD »G JD « E I SD« ECECK1 , EAL , GJ R« EIRL, EC ECK2,! STQP_lfJSTgP2 

COM MON '/RATIO/ A NY R 

DIMENSION DI SS14) »GAMA12( 5 if AL0A0J9) 

NAME L I S T / C Y L D A T A / I S T □ PLiXS TO P^f CA RGE l f CA^C E? f A NX, AN YR t JE L,R> 

_J 1 EX ,EY,GXY, ENUXY,RO 

COM M U N/ LI ST / C ARGE1 ♦ CAR GE2 ; — 

PI=3. 14159265 

2_ FOR MAT I BE 15*6.) _ 

I F ( J J J . N£ . 1 ) _G 0_ TO 1111 

__R E AOI 5 f CYLDAT A_i 

1 1 il__CON li NUE : : ; 

A N Y = ANY R* ANX ; 

E N U Y X= EY* ENUXY/E X 

£NT = j.-ENUXY* ENUYX ' 

QR ( 1 J = E X / ENT 1 - 

QR( 2)=EY/£NT 

j)R( 3 ) -ENU YX»£ X/EN T ____ — 

QR ( 4 ) ~ G X Y __ : 

CALL QL AMN A ( QR ,0» ,QLt4) _ 

CALL QLAMNA(QA,A5_*_iQL t JL2 

CALL ASEMBLI D , QR,A,EAS „ E AR , 32S ♦ B2R, CARGE1 +CARGE2 ) 

; “CAL CULAT I ON OF STRAIN : 

A ( I ) 5 A 1 1 y V EAO 

A { 2 ) = A< 2 ) »EA L 1 

DEL TA=AI1)*A( 2)-A13)**2 

A A 1 = A 1 2 > /DEL T A : 


AA2=A ( 1 J /DELIA _ 

AA? -A ( 3 ) /DE LTA 

AA4 -_ x . / A( 41 

E P S LX= AA1 *AN X + AA3*ANY 

£ PSLY=AA3*ANX*A A 2*ANY 

C * * * * ** * CALCUL ATiON OF STRAIN IN EACH LAYER 

DO 2000 1 = 1,3 

CALL STRAIN! AL PH AXLLi EPSLX » EPSL Y, EP S L It I ) +EPS L2! U ,GAMA12(_IJJ 

_2 000 CONTINUE 

E PS LI !^i = EPSLX 

EPS L2 (4 J = EP SL Y 

G AMA1214I=Q» 

C A L L STRAIN! 45_*_ L E_PS LX,EP SL Yj_E_P SJLH 5JL»JP S L2 J.51_t_G AM A 12J 5JJ„ 

C J 0 F I ND_ L PAD CARRI E D B Y EACH MEMBER ***W*»«*»*« *** 

"DO lo: 1 = 1 , 3 ~ ~ 

__ PXS II ) = EPSLX*EASm 

jo PXR ( 1 )=£P S LY*SAR! I ) 

A ( i ) = A 1 1 ) - EAD 

A<2 )=A<2)-£AL 

P X S K= A ! 1 >*EPSLXfrA (3)* EPSLY 

P Y S K= A ( 3 i > E P S QUA! 2) * E P SLY 

IF! JJ J* EQ# MAX J J J ) PRINT 2 , EA D, G J D , E I SD,E CECK 1 , £AL , G J R , £ I R L , ECECK2 

c ****** TQ FIND BUCKLING L 0 A D $ * ********* * * 

I D0_ 1000 1 = 1 

1 000 D iSSt I J = D ( 1,1) __ 

CAL L I IRA T E_( 1 ,30,0,2 , NC QUNTC 1 ) , BUCKL 1 U , NCR ( 1 ) , HC Rl 1 1 , A, 

i D I SS, EL , ftV ; 

pjl =2.*PI*PI 

IF! i ST 0P2.EQ.1) BUCKL I 2 l = 8UCKL1 1 1 

I F t ISTQP 2 » EQ > 1 ) GO TO 60 00 _ 

E ALjfO •_ ; 

GJR=Q. 

EC E CK.2^0 * ; 

eiRL=0. 

CALL 1TRATE 1 I, 5, 0,2, NCQUNT (2 ) , BUCKL ! 2) ,NCRI 2) «MCR ! 2) , A , D I SSr 

1 ALPHA ( 8 ? ,R ) 

8000 CO NTI NU E 

I F ( I S T QP1 » EQ« 1 ) BUCKL ! 3) — BUCKL 12) 

I F ( I STOP 1 « EQ » 1 ) GO TO 7000 

£ AD=0 » ~ 

gjd1=6. 

ECECKif O. . 

1 Jo I S L)= 0- 


1 800043 
1800044 
1800045 
1800046 
1800047_ 
1800048 
.1300049 
1300050 
1800051 
1800052 
JL 80 0 0 5 3 
1800054 
1800055 
1800056 
1800057 
1800058 
1800059 
1800060 
1800061 
1800062 
1800063 
1800064 
_ie00065_ 
1800066 
J 800067 
1800068 
1800069 
1800070 
1800071 


1800072 

1800073 

1800074 

1800075 


1800076 
1800077 
1800078 
1800079 
_i 800080 
1800081 
J 800082 

_1 800083 
J 800084 
1800085 


N NNN=_INT(P I_*R / ALPHA C 711 

I BUC ( 3 ) gjo'i £2 0 J 

I F ( . *N N . E .Q >_Q„»JD.K>N N NN ♦ EQ . 1 ) G O TO 3 838 

00 .100 _ j= \_i\Q 

NNX=NNNN*J 

_D0 100 I=l,_30 - 

CALL b UC K LG l I .J N :.X_ jAt. OJ S L S_* ANX X »jALJPH Alii: ) ,t„R » 

" I F l AIM X X . L T . 6 U 1C k L( 3) ) MC R (3MJ 

I Ft ANXX • L T 6 UC K L(3)l NC R ( 3 1=_J_ — 

I F ( ANXX * L T.-.BUC K LJ JLI JLJLU C KL iAl-AMX 

100 CONTI NUE 

88 88 I F ( NNNN . EQ.O .OR. NNNN *_EQ »X ) CALL I T R ATE C 1 , 30 , 0 * 2 » NC QUNT ( 3 ) r SUCJ<LJ 

ZI_Il.il NC.R13) ^ rMCR l 3r,At_DiS_S.» AL^A18j_i.R) — 

7000 C ON TI NUE — — 

~ I F 1 ISTQP1 «EQ » 1 ) GO TO 5 000 : : 

1 B UC KL ( 4 ) = P II * ( S Q R r ( D IX? 1 > * 0 ( 1 *_2 U *■ D ( 1 , 3 )>DULt 4.1* ? * U A L P H A JJL 2 1 

BUCKL {5) = PU*(SQRT(D( 2, 1 )*D(2,2) )>D(2,3 ) »2.*P ( 2, 4) ) /B 2$ 

B UC KL ( 6 ) - P I I» t SQR T(Dt3«l)*D13,2 ) i +D (3 ,3 ) »2 . *0 ( 3 , 4 ). 3 / t .8* AL PH A U 2)1 

_5000__CQNTI NUE 

_ l F t I ST0P2 . EQ * 1 ) G O TO 5 001 — 

BUCKL C 7 ) = P1 1* C SQRT ( OU, l)* D(l, 2 ) )>D< 1 , 3 ) »2 . *0 1 1 » 4JJ./ALPHAJ 1J) 

BUCKH 8) =PI I* t SQRTl 0 t4 f 1 ) *0(4,2 ) )>2. *D( 4 r 4^X/8 2R. 

~ BU CKU9) =Pli»( SQRT IOC 5, ll *0(5.2 ii + D 15,3) »2. * D ( 5 , 4 ) ) / I *8 *ALPKA_U7j_l 

50 01 CONTINUE 

__ A LOAD ( 1 ) = ANX 

ALO AD l 2 )- ANX_ 

A LOAD ( 3 ) = P X S K : 

A LOAD (4 ) * P XSK - 

A LO AD ( 5 ) = PX $ ( 2 ) 

ALO ADI 6 ) = P X S { 3J 

__ ALO AD i 7) ~PYS K Ll 

ALO AOiaJrPXR.t 2i 

ALO AD ( 9 ) -=PXR ( 3J 

C MAT E ft I At PROPERTIE S USED ARE FOR THORNELL 300 NARMCO 5208 

C : — 

C EXX 2.12E 7, E22 2«39£6 , G12 6.5E5, NUE12 » 31 . 


DO 33 00 lz_l_*_9 

FUNCTNt Ii=SlBAR(U : _ 

I F f ALOAOt I>»LT»0«i GO TO 3300 

iVt A LO AD ( I ) > GT » BUCKH I j JFUNCTNt I) * ( ALO AD ( U - BUCKL U ) )/ AL OAD ( I ) 

3 3 00_ _C 0 N T I N U E — - 

DO 3400 1=1 0,14 


1 yuuuyb 


180008 7. 
_1 800088 
_1 800089 
800090 
_1 800091 
_1 800092 
.1800093 
_ 1800094 _ 
1800095. 
_1 800096 
1800097 
_180009.8 
1800099 
1800100 
1800101 
1800102 


1800103 
180 01 04 
1800 105 
1800 1 06 
_1 800 1 07 
_1 800 108 
_1800 1 09 
_ 1 800 110 
_ 1 800 111 
_ 1 800112 
1800113 
_ 1 80 0 1 1 4 
_1800 115 
_1 8001 16_ 
__1 8001 1 7_ 
1800118 
1800119 

__i 8001 20 
1 800121 
_1 800122 
1800 123^ 
__i 800124 
800125 

18001 26_ 

1800 1 27 

1800128 


FUNCTN!I)=$IBAR( i) 

1800129 


FUN' T N(U5) = SIBAR (1*5) 

1800130 


FUN-«Nl I*10)=SIBAR( I *10) 

1800131 


IF! EP $L 1 ( 1-9) .LT.O.JGO TO 3420 

1300132' 


IF ! ABSiEPSLll I-9H .GT. .012 ) FUNCTNi 1 ) =IEPSL1 ( 1-9,)-. 012 j /EPS! 1 1 1-9 i 

1800133 


IF ( ABS! EPSL2 (1-9 M • GT • • 02 ) FUNCTN ( 1*5 ) = ! EPSL 2 ( l- 9 )— .02 ) /EP SL2( 1-9) 

1800134 


I F l ABS ( GAMA1 2 l 1-9 )) . GT . .01 5) FUNCTN! I* 10) = (G AM A1 2 ( I-9>- .0 15 ) /GAMA12 

1800135 


1(1-9) 

1800136 


GO TO 3400 

1800137 


3420 I F ( ABS ! EPSL1 ( 1-9) ) .GT. .01 )FUNCTN ( I ) = { EPSL1 < 1-9 )- .0 1 ) /EPSL 1 U-9) 

1800138 


I F( ABS l EP SL 2 U -9)) .GT. • 0045) FUNC TN( I>5) = ! EPSL2( 1-9)-. 0045 ) /EPSL2 ( I 

1800139 


1-9) 

1800140 


IFt ABSIGAMA12U-9) ). GT . .0 1 5) FUNCTN! I> 10) * (GAMA12 ( 1-9) ~. 0 15 ) /GAMAi 2 

1800141 


1(1-9) 

1800142 


3400 CONTINUE 

1800143 


FUNCTN125)=2.*PI*R*EL*(TS+ASS/ALPHA(7)+A$R/ALPHA!8) ) 

1800144 


IFt JJJ.NE.MAXJJJ) GO TO 1113 

1800145 


8 FORMAT!//********* SOME RESULTS *******///) 

1800146 


PRINT Q 

1800147 


3 FORMAT ( * NX NY/NX EL Rfr) 

1800148 


PRINT 3 

1800149 


PRINT 2 ,ANX # ANYR, EL » R 

1800150 


4 FORMAT ( /* EX EY GXY ENUXY*) 

1800151 


PRI NT 4 

1800152 


PRINT 2, EX » EY * GXY * ENOXY 

1800153 


5 FORMAT ! //* .ALPHA1 ALPHA2 ALPHA3 T1 

1800154 


1 T 2 T3 ELS * ELR*) 

1800155 


PRINT 5 

1800156 


PRI NT 2 , ( ALPHA ( I) , 1*1,8) 

1800157 


6 FORMAT!/* T1S T2S T3S BS 

1800158 


1 HS*) 

1800159 


PRINT 6 

1800160 


PRINT 2,( ALPHA!!) ,1*9,13) 

1800161 


7 FORMAT!/* T 1R T2R T3R BR 

1800162 


1 HR*) 

1800163 


PRINT 7 

1800164 


PRINT 2. (ALPHA! I) #1=14,18) 

1800165 


2001 FORMAT!//* STRAININ THE FIBER DIRECTION *) . 

1800166 


PRINT 2001 

1800167 


2003 FORMAT!* ALPHA EPSL1 EPSL2 GAMA 12 ,*) 

1800168 


PRINT 2003 

1800169 


2002 FORMAT !F7. 3 , 4X, 3 El 5. 6/1 

1800170 


ZE=0 

1800171 



ZEE~45. 



1800172 . 


PRI ' 2002 ? l ( ALPHAl I ) * EPSLli I) » EPSL2! I ) ,GAMA1 1)1. 1 = 1. 3) 



1800173 


- 2 go 2 . Z £ t E PS L 1 1 4 1 ,EP SL2 (4) * GAM A 1 2 ( 4 ) 



1800174 


dot wt ?nn? - 7 P F . FP<;i 1 I S ) „ FPSl ?!5).GAMA12!5) 



1800175 


1123 FORMAT (* STIFFNER PROP ERT I ES*//*EAD=* ,E 14.7j *GJ D=* * E14 .7.!*E I SD=* 



1800 176 


1 f E14.7,*ECECKl=*, E14.7) - - 



1800177 


PRINT 1 123* EAD , CJL) . E I SD *ECECK1 



1800178 


9 FORMAT!/* 1 AXIAL STRAIN TRANSVERSE STRAIN*) . .... _ 



1800179 


PRINT 9 



1800180 


PRINT 2 .EPS! X r FPSl Y 



1800181 


20 FORMAT l***LGAD CARRIED 3Y EACH MEMBER*) 



1800182 


PRI NT 20 



1800183 


11 FORMAT ( * PXS.PXR.PXSK.PYSK*) 



1800184 


PRINT 11 



18001 85 


PRINT 2, (PXS( I ) tl=lt 3) 



1800186 


PRINT 2 * ( PXR ( I ) *1-1.3) 



1800187 


PRINT 2.PXSK.PYSK 



1800188 


12 FORMAT! //*♦*** BUCKLING LOADS*) 



1800189 


PRINT 12 



1800190 


13 FORMAT 1 /*GRQSS BUCKLING PANEL BUCKLING SKIN BUCKLING*) 



180019 1 


PRINT 13 



1800192 


PRINT 2. ( 8UCKL I I ) *1 = 1*3) .. ~ - 



1800193 


K ~FORMAT(/*LOCAL BUCKLING STIFFNER*! ... 



1800194 


PRINT 14 



1800195 

..... 

PRINT 2 »!BUCKL( I) * 1=4.6) 



1800196 


15 FORMAT!/* LOCAL BUCKLING RING*) 



1800197 


PRINT 15 



1800198 


PRINT 2* < BUCKL! I) • 1=7, 9) _ .... 



1800199 


16 FORMAT!/*** M CR ( l i * NCR ( 1) * NC OUNT II ) *) 



1800200 


PRINT 16 - - 



1800201 


21 FORMAT (31 8) 



1800202 



PRINT 21. ( IMCR! I) fNCR(I) .NCOUNTU) ), 1-1*3) 



1800203 


CALL OUTPUT I FUNCTN tALPKA» 1 ST0P1 » IST0P2 «EL«R» RO* ANX ) 



1800204 


1113 CONTINUE 


— 

1800205 

1800206 


RETURN .... — 

ENO ... — 



1800207 



N E A SFMBMp* Q R»A»E AS»£ARyB 2 St B 2 Rf C A R G 
COM . .jN /AREA J f S , A_SS_f ASR 

“DIM e ns ion 
1 ,A 14 L 


,CARGE2J 


q r ( <V ) ♦ THIC KJ 6 ) « KP ( 6 ) » PS t 4 ) « D( 5 * 4 ) » DSSL2jrAV»^ Lg_EARj;3 j 


COMMON / Q / Q L ( 4 * 5 J - ■ — — ^ — ~ t c r n oT 

COMMON/ S A VE1 /EADf GJ 0 «E ISO » ECECK1 » E A L » G JR» E IRLyEC E CK2»I STOP! f I ST OP 2_ 

COMMON /AESOPD/ A DA T A< 5QQQ_) 

REAL AL P HA l 2 0 J J _ 

EQUIVALENCE! A DAT A ( 742 ) * ALPHA) 

00 1 i = if 3 


CALL QLAMNAlQRtAL PH A t l ) t QL* 1) 

"T HICK l I > = ALPHA V_l f3J 

THI CK l 7-1)= ALP HA! I »3 ) 

KPCiR " 

KP! I*3)=4-J 

CONTINUE 


pis T=A L PH At 4 ) + ALPHA (5) + AL P H A ( 6 ) 
CAL L 
DO 2 

'oil \ i) = osai 

TS=2.j*DI ST_ 

G X Y $k = A ( 4 > / T S 


ST I FF ( A, K P YoS t T H I CKt 0 IS T» 6 i 

I - 1 1 4 


l F t I S TOP X . EQ . 1 ) GO TO 50 00 


CA L L ST I F PR ( AL PH A ( 9 )»ALPH A(10)»ALPHA( l l)fALPH At 12 ItALPHAtiaJ » T S ? 


E A D » G J Q, EC E C K 1 ,E I SP » DSS * ALPHA! 7 ) i E AS * 82 S »C ARG.El _y_ASSf.GXYS Kt GXYBSTJ^ 

DO 3 1 = 1 1 4 — 

D ! 2 ? I > - 0 S S l 1 1 1 ) 


3 D( 3 1 1) - DS S! 2 t I ) 


XI YF AKThTT* AL ? HAt 12 ) *GX Y3 ST/ALPH A 1 7) 


5000 CONTINUE _ 

IP( IST0P2*EQ.I) GO TO 5001 

r ALL STIF PRTaLP HA ( 1 4 )T A L P HA ( 15 ) » AL PH A 1 16 ) t ALP HA (1 7 ) > ALPHA ( 18) i 


T S VEAL', G j R7E C £ C< 2 1 E 1 R L > D S S, A L PH A 1 8 ) »£ AR T B2R » CARGE2 » ASRt GXYS KtjjO^ 
BR N ) — 


DO 4 1=1*4 


ot 4 , n = Dssu»n 


4 D 1 5 tT) =DSS t 2 * 1 ) 


A'C4) = AlV)' + l7T^LPHAt 17~)*GXYSRN/ ALPHA! 8) 


5001" CONTINUE 
RETURN 
END : 


1800208 
1800 209 
1800210 
1800211 
1800212 
1800213 
1800214 
1800215 
1 800216 
1800217 
1800218 
“1800219 
"1800220 
__1 600221 
1800222 
1800223 


1800224 
1800225 
”1800226 
1800227 
1600228 
1800229 
1800230 
1800231 
1800 2 3 2““ 
1 800233 
1800234 
J. 800235 
1800236 
~1 800237 
_1 800238 
_1 800239 
1800240 
JL800 241 
_1 800242 
_1 800243 
1800244 


1800245 
1800246 
1800247 
18002 48 


BLC ' DATA 

1800249 


DIM I ON BUCKL(9)*EAS(3J »EAR(3) 

1800250 


COMMON /DATA/ BUCKL , E AS , E AR , B2S , B2R 

1800251 


CaMMON/SAVEl/EAD,GJD,EISD»ECECKUEALtGJR,EIRL*ECECK2, ISTGP1, 1STOP2 

1800252 


COMMON /AREA/ TS»A$S,ASR 

1800253 


COMMON /MAT/ EX , EY , ENUXY , GXY , RO 

1800254 


COM MON/ LIST/ C ARGE1 » CARGE2 

1800255 


DATA IS TOPI, I$TGP2/2*0/ ,CARGE1, CARGE2/2* W, EX/10. E6/ , E Y/10« E6/ , 

1800256 


1 GXY/3. 75E6/, ENUXY/. 33333333 /,RO/.l/ 

1800257 


DATA TS,A$S,A$R/3*0./ 

1800258 


DATA BUCKL /9* 1 *E 307 , E A S/ 3*0 . / , E AR/3*0. / , B2 S * B2R , EAD, G J 0, E ISD, ECEC 

1800259 


i K1 ,EAL,GJR,EIRL,ECECK2 /10*0.A 

1800260 


END 

1800261 



SUB' ^UTIN E BUCKLGlM,N,A,D,ANX,Et.,FU . ' 1800262 



COM..jN/SAVEI/EAD f GJD# E i SD t ECECK1 ,EALf GJR# El RL # cCECK2» I 5T0P1 * I STOP 2 

1800263 



COMMON /RATIO/ AN YR 

1800264 



DIMENSION A(4) f D(4) 

1800265 



P I® 3# 141 5927 

1800266 



PIM=FLOAT (MJ#PI/EL 

1800267 



RN= FLOAT ( N> /R 

1800268 



5 A 1 1 = t At 1 J + E AD ) *P IM**2 + A.( 4 ) *RN**2 

1800269 



A12 = l A i 3 ) *A 14 ) )*PIM*RN " " 

1800270 



A13-A ( 3 ) *P I M/R +■ EC ECKi *P I M**3 

1800271 



A22 =A(4) *PIM**2M A(2 J +EAL } *RN**2 

1800272 



A23®CAl2|fEAL)*RN/RfECECk2*RN**3 

1800273 


• 

A 3 3 = ( Dt 1 ) 4-EISq ) *P IM**4^ (2.«t D( 3 )F2.*D{ 4) ) *G JD+G JR ) *PI M**2*RN**2 

1800274 



1M 0 (2 ) + E I RL j * R N*4~4 +• l A t 2 ) + E AL J/R**2 *2. *ECECK2*RN**2/R 

1800275 



_APP®All*A2 2-A12*A12 

1800276 



APX = A33M A12*A23-A13*A22)*A13/APPM A12*A13~All*A23i*A23/APP 

1800277 



ANX®APX7<PIM**2*ANYR*RN**21 

1800278 



RETURN 

1800279 



1800280 


SUB " I UT1 NE XTRATE( MIfMF»NI t ND£L f NCOUNT f C M IN y N' ,MCR «A 

"DIMw.4 s i ON “a ( 4) »pi 4) I : 

F OR MAT ( I HO »» WARNIN G VALUE OF N WENT BE YON D 100 g7_/ j 

NCUUNT^Q ~~~ 

NN=NI _ 

NDE LT A= NOEL 

D 0 10 0 1 = M I » MJF 

I F ( l • N E • 1 ) " NN=NC 1 __ 

IF ( I.NE.l) NDELT A= 1 

INN = L : 

_NCH ECK=0 

i F 1 NN .GT. 100 ) P RINT 1 

I F { NN.GT.lo6) GO TO 100 

CALL BUCK LG II , NN , A, D , FUNCTN, EL » RJ 

“NCOUNT = NCOUNT + 1 

JFt INN. E Q. 1 ) GO TO 30 

IF l F UN C T N . L T . C F U NIN) GO TO 30 

NC H E C K = N C H E C K Vf . 

' NOELTA=-i 

'NN=NN*NUE_LTA 

INN =2- 2_I 

IF! NCHECK. E'Q . 2 ) NC 1= NN+2 . 

IF( NCHECK.EQ .21GQ T 3 IQ 

C FUN T N=FU NCTN ~ 

GO TO 50 

1 F ( I. N E.MI ) GO TO 200 

3 c r =m r __ 

NCR=NC 1 

I F { I.EQ.MI) GO TO 70 

IF l CM I N.LT . CFUNTN ) C~FU NT; N=CM I N 

■ i F ICFUNT N • L f ."CM I N ) M CR=1~ ' 

I F { CFUNTN.LT* CM IN) NCR = NC1 

_CM I N=CFU NTN 

CONT INUE 

' RETURN 

E ND 


? MCR* A » Of EL »R) 



JSU8 ' ^UTIN E QLAMNA1 A t THETA»3fK.) 

Dlh^./SION A (4 > y 8t 4 > 5 1 ; 

T HE TE= THE T A*3. 141 59265/ 1 80 . 

JS=SJ N { THETE) 

C=COS( THETE) __________ 

C 4=s"c** 4 

S4=S**4 

C S2 2= C *C* S*S 

C$3=C*S**~3 

_S C 2_f s* C* *_3 ; 

Ai = 2 • * ( A l 3 ) fr2«*A (4H*C S22 

A2= A U ) - A 1 3 )-2.*Ai_4j 

A3= A2 *- A ( 2 J-A13J 

"A4= At 3>-Al2)J-2.*A(j4) 

B ( 1 yK)=Am»C4 + Aj>A( 2)* S4 

BJ2 ,K) = A( l)» S4f Ai»A( 2>*C4 

B(3,K) _= ( A ( 1 ) f A( 2) -4.*A(4 )]*CS22«-Am *(C4+S4) 

B 14,K ) = A 3 * C S 2 2 i- A t 4 ) * ( C 4 <- S 4 ) 

RETURN 

~END " " __ 


1 8003 17 

1 80.031 8_ 

1800319 

1800320_ 

1800321 

_1800322_ 

1800323 

1800324 

1800325 

1800326. 

1800327 

1800328 

1800329 

1800330 

1600331 

1800 3 32_ 

1800333 

1800334 

1800335 

1800336 


SUB WINE STIFF! A, <P,D,T H ICK, 0 1 ST, NU 

COM..JN /Q/ QU4,_51 

blM E NS I ON A (4 J , < P < 6 ) ,0 ( 4 j > TH ICK! 6) 

DO 50 l = i' , 4 

A ( [ ) = 0 ~. 0 

b ( i *=o.o 

50 CONTINUE 

HK1=-D I ST 

DO" 100 J*itNL ; 

KK=_KP_l U 

HK2=HKl*THlCta_U _ 

HP A=HK2“HK 1 

"" HPD = (HK2**3-HK 1**3)/3. 

__ DO 20 K= 1 >4 

A(KJ=A(K)*QL (K, KK)*HP A 

20 U (K)^D( K)»Ql(K f KK)*HPD 

HkT*HKF __ 

10 OT C 0 NTINUE " 

RETURN 

EN D I 


1800337 
_1 80033 8_ 
1800339 
18003 40 
1800341 
1800342 
1800343 
1800344 
1800345 
1800346^ 
“l 80034 7 
1800348 
1800349 
1800 350 
1800351 
_1800352__ 
Jl 800353. 
_1 800 3 54 
_1 800355“ 
1800356 


SUB 'JUNE STIFPFUT1,T2,T3,BB,H,T.S,£AD,GJD,EC' U , EISD, DS , ELD , EA, B 

1800357 


1 AB, v..iARGE,ASSS,GXYSK,GXYB J 

1800358 


DIMENSION THICK I 6) , KP < 6 ) , A( 4} , Dt 4 ) , EX (3 J *GXU3) ,DSL2 . 4 > , EA 1 3 ) , B 13 } 

1 80 0 3 5 9 


COMMON /Q/ QL(4,5J 

1800360 


THICK(1>=T2 

1800361 


kp ( n=5 

1800362 


0 l S T= T2/2* 

1800363 


CALL STIFF<A,KP,D, THICK, DIST,1> 

1300364 


DO 1 1 = 1,4 . 

1800365 


1 DS( 1, li=OIIJ 

1800366 


EX 12> = 1 Al 1)*A(21-A( 3)**2) / i A 1 2)*T2J 

1800367 


GXY 12)=A(4)/T2 

1800368 


THlCKt l)=T2/2. 

.1800 369 


THICK(2)=T3 

1800370 


THlCK(3)-=T2/2* 

1800371 


01 S T= [ T3 + T2 ) /Z* 

1800372 


K P ( 1 > = 5 

1800373 


KP( 2)=4 

1800374 


KP ( 3 ) = 5 

1800375 


CALL STIFfM A , KP , 0 , TH l CK , D I ST ,3 ) 

1800376 


DO 2 i=l»4 

1800377 


2 DSl 2,1) =0(1) 

1800378 


EX(3) = UU)*A(2)-A(3>#*2)/IA(2)*(T2 + T3* i 

1800379 


GXY(3)=A(4)/{ T 2 ♦ T 3 ) 

1800380 


EZ= (Ql(l,4)*QL (2,4)-Ql(3,4)**2)/QL(2,4) 

1800361 


EF= (QU 1, 5>*QL( 2»5)-QLC3t 5)#*2J /GU2,5i 

1800382 


EX (!> = ( EZ*T1*EF*T 21/ (T1+T2) 

1800383 


GXY (l)MQL(4,4)*Tl*QL(4f5) J *T2)/ 411*121 

1800384 


BETA=ATAN(.1*8B/H) 

1800385 


TAVG*t fS*Tl*T2)/2. 

1800386 


B(18=.4*BB 

1800387 


B i 2 )= (H-TAVG) /COS( BETA) 

1800388 


B i 3 ) = .8*BB 

1800389 


A~SSS=BU)*IT1*T2)*2.+2.*B<2)*T2*B(3)*(T2 + T3) 

1800390 


BAB=B(2* 

1800391 


EA(l) = EX(l)*8(l)*(mT2> 

1800392 


EA12 )=EX(21*B(2)*T2 

1800393 


EA( 3J=EX<3)*6(3)*<T2*T3> 

1800394 


E AA =2 • * E A ( 1 ) + E A ( 2 ) *2.+EA(31 

1800395 


ZBAR=(2.*EA( 1 )*TAVG + EAt 2) *IH*TAVG) *EA13) *HI/EAA _ 

1800396 


Z8AR=ZBAR*CHARGE 

1800397 


E IC = (EXC1 J*Bl 1)*IT1«-T2 i ) / 6 . + I EX { 2 > *T2 * ( 8 ( 2)*C05(8ETA ) )**3 

1800398 



)/ 16 • *COS ( BET A ) )_M EX( 3J *B J 3 1 * 1 T 1»T 

“V 2 i A 12) * U HtlAV 

GJ2=oiZ)/ (I2*GXYjijJ. 

G J 3-B l 3 ) / { l T2_4-T3 )j?_GX Y.L3JJ — 

~G 14-B B/ (TS^GXYSK) — 

GJ4=0. 


+2-*EA ( 1 )»T AYG** 2 


G J = ( 4 . *J * 9 *BB* H )* * 2) / ( 2 . * G J2 «- G J 3 »G J 4 ) * t L1±T 2 1» 

G X YB^GXYl?)* T 2 * ( ♦ 3*BB»H»T AN( BET A1 ) / ( » 8 *B Bjji /C^SXB.ElAj_) 

E AO = E A A / E LD 

_GJD =GJ/ EL D ; 1 

E I S D= E I C/ EL D _ 


ECECKl = eAD?_Z.BAfL 

RETURN 

END 


.1800399 
_1 800400.. 
.1 800401 
1800402 
JL 800 40.3- 
_1 800404. 
1800405.. 
1800 406. 
..1800407. 
..1800408. 
.1 30 04Q9_ 
_1 80041 0 
.1800411 
_180Q412 
1800413 



SUB IT I NF ST RAIN 1 T H £ T A , £ PSULtJ^S 

COMhuN /AES DP D / AOATAC 5 ,0.001 

RE aL.AL.PH A (20) — 


_EQU IV AL ENCE ( _ADA1A (, 7_4_2 1 > ALP. HO. 

f HE TE = 3 . 1 41 5 92 6 5* THET A/1 80* 

OCOS(THETE) 

""$* S INC TH ET E i _ 

CC=C*C 


EPSLl=CC^EP_Sl x» ss*e ps.ly 

E P $ L2j= S_S* EPS LJL±CC*JP SJJL 

GAM A1 2=2 . *( ~CS*EP S LX +C S*EPSLY_1 

return” 

END 


%GAMA12) 1800414. 

180041 5_ 

180041 6_ 

~ “ 1800 41 7_ 

“““ " 1800418 

" 1800419. 

— 1300420. 

1 800 42. 1_ 

~ ' 1800 422 

“ 1 80042 3_ 

“ “ 1 800424. 

“ “ 1 800425- 

“ 1800 426. 

~ 1 80042 7.. 

~~ 1800428. 


SUB' VTIN6 OUTPUT { FUNCTNf ALPHA t l STOP! t 1STQP2 ,/ R>ROt ANX | 1800429. 

C THIS. SUBROUTINE OUT PUT S THE OPT IMi.ZAl.ON RE SUITS _ 1800430 _ 

___ j) l MEN SI ON FU NCT N( IOO) » AL PH A_( 1 00 J 1800 43 l_ 

1 FORMAT 1 1 HI 1800432„ 

Z FORMAT ( * OPTIMIZATION RESULTS FOR UNSTIFFENED CYLINDER*) 1800433 

3 FORMAT I *" OPT I MI ZAT ION RESULT S F OR STIFFEN CYLINDER WITH NO R ING*) 18QQ434 

4 FORMAT I * OPTIMIZATION RESULT S..FOR STIFFENED CYLINDERS) 1800435., 

3_ FORMAT (* BUCKLING CONSTRAINT. VO ILATIQNS, GROSS., P ANEL, SKIN BUCKL* ) I800436_ 

6 FORMAT I * LOCAL BUCKLING C GNST RA I NT_ VO XL AT I ONS OF ST IFFNERS*) 1800437.. 

7 F ORMAT (» LOCA L BUCKLING CGNSTRAINf VO I L AT IONS OF RINGS*) 1 80043 8_. 

_8_ FORMAT i 1 H 0 ) _ 180043 9 _ 

9 FORMAT! / /IX > 5F22*6 / ) 1800440., 

10 FOR MAT l * CONSTRAIN VOILAT IONS FOR STRAINS, E P SL 1 < I) » EP5L2 ( 13 >GAMA1 1800441, 

12 ( II*) ™HT„™ " - 1800442,_ 

ii FORMAT (V jBAR7R = *tF8 3t5_X f *NX/R= * , F 8 . 3 , 5X **TB ARD/R=* , E I 5. 8 1 1800 443_ 

~\Z FOR MAT l *_ALPHA 1 =*F 5. 1 , 3 X> * AL FHA 2=* ,F 5 , 1 , 3X, » ALPHA3 = * , F5 * 1 , 3X ,* T l~ 1800444, 

1* ^F6 *4 ^ 3X. » *T2 =7 fF6*4»3Xf * 1800445., 

13 FORMAT ( * T 1 S = * , F 6 . 4 , 3 X , * T 2 S_= * , F_6 . _4 f 3 X ,*T3 S= * , F 6 . 4 > 3X, * B$- » F5 >2 t 3X» _ 1800 446, 

2 2*HS ss *tF5*2 r5X»*ELS=*» F6 » 2 3 ; 1 800447,. 

14 FORMAT ( * T 1 R = * , F 6 • 4 »_3X » *T2R=*,F 6 .4 f 3 X , * T3 R= * » F 6.4* 3 X t AB 5 .* 2 , 3 X 1800448 , 

~~ 2**HR=*,F5*2 r 5X**ELR=*,F6.2) 1800449, 

15 F OK MA T ( * NX^*,F8. 3» 5 X , *L = *t~F8.3 t 5X » *R=* >F8» 3 > 1 8 0 0 4 5 0_ 

PRINT 1 , 1800451... 

IF( IST0P1. AND. isI0P2.EQ,lJ PRINT 2 _ _ 1800452. 

I F ( IST0P1 * E Q • 0 I PRINJ __4 ‘ 1800453, 

__ I F { ISTOP2.EQ.O)PRI~NT _3 1800454 

J PfU NT 8 1800455_ 

PRI NT 1 5* ANX t EL >R 1 800456, 

P RI NT 8 1 800 45 7 .. 

__ PRI NT 5 1800458, 
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APPENDIX B 


MODIFIED STIFFNESS MATRICES 


FOR USE OF BUCLAP 2 




SECTION B 


EFFECTS OF BOUNDARY CONDITIONS 
ON THE BUCKLING OF AXIALLY 
COMPRESSED CYLINDRICAL SHELLS* 


This section will form the basis for a paper to be submitted for pub- 


lication in the AIAA Journal. 




BUCLASP 2 requires the input of [A], [B] and [D] matrices (see Ref. 6 
for definitions) which , for the present problem y are given by 




1, •' INTRODUCTION . . 

1.1 PRELIMINARY CONSIDERATIONS 

The literature is replete with investigations devoted to 
the buckling of unstiffened, isotropic, complete cylindrical shells 
under axial compression. In contrast, relatively little attention 
has been focused on the corresponding buckling problem for cylin- 
drical panels.'. However, unstiffened, isotropic cylindrical panels 
are frequently employed in a vast number of structures, such as, 
for example, in launch vehicles. Furthermore, knowledge of panel 
buckling loads is needed for consideration of the local panel 
buckling modes in the analysis and minimum-weight design of stringer 
stiffened cylindrical shells. Therefore, the present paper is de- 
voted to a study of the buckling behavior of unstiffened, elastic, 
isotropic, cylindrical panels. The loading condition of uniform 
axial compression is chosen as this condition is often the critical 
one, especially in aerospace applications. Buckling loads are 
presented for panels with eight sets of boundary conditions along 
the straight edges of the panel. Four, sets are considered for both 
simply supported (w ~ M = 0) and clamped (w - w, - 0) straight 

J ■* 

edges* The eight sets of boundary conditions are designated by 
SSI, SS4, CC1, CC 4 and are defined below. The boundary 

conditions for the simply supported straight edges are 


SSI 

: w 

= M 

= N 

II 
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ii 

o 
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xy 



SS2 

: w 

= My 

= N xy 

= v ■ = 0 
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(classical) SS3 

: w 

= M 

y 

= u 

= N =0 
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SS4 
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= M 
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= u 

ii 
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ii 
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The corresponding boundary conditions for the clamped straight 

edges are , 

CC1 
CC2 


CC3 

CC4 


V 7 = v/ — N — N — 0 

,y *y y 

w = w = N = v = 0 

w ,y xy 

w = w ,y = u = N y = 0 

w = w = u =V=0 


( 2 ) 


In these equations, u, v, w and N y , N xy , M y are, respectively, 
incremental displacements and stress resultants that take place 
during buckling. Their positve directions are shown in Figs. 1 and 
2 along with the other stress resultants considered herein. Different 
sets of boundary conditions along the curved edges of the panel will 
not be considered here since it seems reasonable to conjecture, as 
was done in Ref. 1, that the qualitative effects of these boundary 
conditions should not differ appreciably from those found for complete 
cylinders, which have been thoroughly documented in the literature. 
Therefore , only one set of boundary conditions for the curved edges 
is considered herein; namely ••classical" simple support edge con- 
ditions defined by 


w 


= M /+>N ~ V = 0 




(3) 




1. 2 PREVIOUS INVESTIGATIONS 

Analyses of the axial buckling behavior of cylindrical panels 

have been performed in a number of references. The references cited 
herein pertain to long and narrow panels analyzed by Donnell-type 
shell theory. In all cases, the effects of prebuckling deformations 


2 



are omitted and the only non-zero prebuckling stress resultant is 


the axial one, N , which is assumed to be constant. The effects 
of initial imperfections are not considered. A. brief discussion of 
the results of some of the other investigations follows. 

Marguerre (Ref. 2) presented solutions for panels with the 
following types of edge conditions: SS2, SS3, SS4, CC2, and CC3. 

The solutions are approximate except for the case of classical simple 
support conditions, SS3, for which the following closed form solution 
was obtained (also presented in Refs. 3-5): 

2 <V K 

<V , n + f ) - <v + | WT ^ 1 ' K .< 2 

pi pi x 


N 


Pi 


XO 


(4) 


(N ) 

x cyl 


K > 2 


In this equation, N is the panel buckling load; (N_) is the 

xo .pi 

classical uniaxial compressive buckling load for a long flat plate 

with simple support conditions on the long unloaded edges, i.e*. 


(N x > 


pl 


. . 2„ . ,3 

tt E ,h x 

9 ' 2 ^ 

3(l-v ) 


(5) 


where h and b are the plate thickness and width, respectively, E is 

Young 1 s modulus, and v is Poisson* s ratio? (N ) is the classical 

cyl 

axial compressive buckling load for a long complete cylinder, i.e.. 


(N ) 
x 


cyl 


(3 (l-v 2 )! 1 ^ 2 


h 2 


( 6 ) 
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where a is the radius of the cylinder and h is its thickness; K 
is a panel curvature parameter defined by 


K = 


(N ) o 1/2 

x cyl _ 1 [3 (1-v ) ] 

TO “I 


x 


,2 

(— ) 
v ah ; 


Pi 


(7) 


9 1/2 

_ 1 [3 (1-v ) 3 2 

- T ( E° ♦o 

TT 


where (see Fig. 1) a is the radius of the panel, h is the panel 

thickness, <j> is the central angle of the panel, and b = a<f> Q is the 

panel width. With regard to Eq. (4) it is noted that Sullins, Smith 

and Spier (Ref. 6) , employed this equation, referred to therein as 

the "Schapitz criterion", in conjunction with a "knockdown" factor 

(from Ref. 7) for (N } and obtained a design curve which provided 

X cyl 

a rather close lower bound to experimental results for unstiffened 
isotropic cylindrical panels. 

For all sets of boundary conditions, Marguerre's numerical 

results show that the panel buckling loads tend to monotonically 

approach a lower bound asymptote, which is the complete cylinder 

buckling load, as the curvature parameter K is increased. For the 

range of K values considered, namely K _< 9 , the complete cylinder 

buckling load was actually reached only for the classical simple 

support set of boundary conditions, SS3 ? for which N x = (N x ) for 

cyl 

K > 2, as shown by Eq. (4) . 

Rehfield and Hallauer (Ref. 8) presented buckling loads appro- 
priate to the eight sets of boundary conditions defined by Eqs . (1) 

and (2). Thus, in addition to the panels studied by Marguerre, Ref. 
8, considered SSl, CC1 and CC4 panels. With the assumptions of a 


RECEDING PAGE BLANK 


not filmed 
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linear membrane prebuckling solution, and sinusoidal axial buckling 

modes, the partial differential ^governing the buckling problem were 

A 

reduced to a set of ordinary linear differential equations with 

constant coefficients and "exact" results were obtained in a straight 

forward manner through an iterative numerical solution of a tran- 

scendential eigenvalue equation. The solutions of P.ef. 8 are based 

on the assumption that N ' > (N ) . This assumption was imposed 

xo cyl 

to insure that the roots of a characteristic equation would always 
be real or purely imaginary, but never complex. As a consequence 
of this assumption, the results of Bef. 8 are somewhat incomplete 
for certain cases, as will be discussed in what follows. 

* 

The results presented by Behfield and Hallauer are replotted 
in Fig. 3 using a different ordinate, namely the non-dimensional load 
parameter . 





N 


xo 


(N ) 
x 


cyl 


( 8 ) , 


Note that buckling loads for the cases SS2 and CC2 are not shown since, 
to the scale of Fig. 3, they are always close to the buckling loads 
for SS4 and CC4 panels, respectively. For SS2 , SS4 , CC2, CC3 panels, 
Marquerre*s (Ref. 2) approximate results are in close agreement with 
the exact results obtained by Behfield and Hallauer. For SS3 panels the 
results of 'Refs'.' 2 and 8 are identical. 

Bef. 8 plots N xo /(N x ) against K and shows complete curves for 
the range of K from K = 0 to K ■ 5. 
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For SS3 panels, it is seen from Fig. 3 that the complete 
cylinder buckling load is achieved for a relatively small value 
of the curvature parameter, namely K = 2 (see also Eq. (4) ) , at 
which the slope of the p - K curve is horizontal. This is hardly 
surprising since the SS3 boundary conditions (w - M y = = u = 0) 

are precisely the conditions satisfied on axial nodal lines of the 
nonsymmetric buckle pattern of the complete cylinder. Thus, for 
K >• 2, the panel is sufficiently wide to permit the formulation 
of one of the (infinite number of) buckle patterns that are possi- 
ble for the complete cylinder at the classical cylinder load (Eq. 

(6)). It is also interesting to note that it can be shown that 
K = 2 corresponds to a panel width b equal to a full wave length 


appropriate to the axi-symmetric buckling mode of the complete 

_ _ , 

cylinder at (see Ref. 3). For CCl panels, p - 1 was 

cyl . _ 

reached within the small range of K values considered in Ref* o. 

Hallauer (Ref. 5) give the following closed form solution for CCl 
panels : 


P = 


. B 

C 

K 


' K 
4B 


K < 2B 


(9) 


where 


B = 1.7428 


( 10 ) 


is the uniaxial buckling coefficient for infinitely long 
with clamped unloaded long edges (Ref. 3). For K > 2B^, 


flat plates 
Hallauer? 
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on the one hand states that p = 1/ and, on the other hand, points 
out that in view of some low buckling loads obtained for SSI panels 
(see Pope’s results in Fig. 3), that there is a possibility that 
the CC1 panels have solutions p < 1 for K > 2B w . However, Eq. 

(9) shows that the P ccl ~ K curve (see also Fig. 3) has a horizontal 
tangent at K = 21^ which strongly suggest that p ccl = 1 for all 
K > 233 (as will, be demonstrated in the present paper). For all 

co 

other cases (with the exception of SSI), the curves in Fig. 3 tend 
to suggest that p 1 asymptotically with increasing K. Now ono 
would expect intuitively, perhaps, that the buckling load for a 
sufficiently wide panel (<j> o <_ 2 tt) , with appropriate support con- 
ditions along the straight edges, should not differ, appreciably 
from the buckling load of the corresponding complete cylinder. The 
results of Kef. 8, since they are limited to rather narrow panels. 

(K < 5)* cannot predict in all cases the value of <f> 0 above which 
the panel and cylinder buckling loads coincide. The behavior of 
p with increasing and arbitrary K will be studied in the present 

work. 

Figure 3 also shows the results obtained by Rehfield and 
Hallauer for SSI panels. As may be seen from the figure, their 
SSI results are restricted to the very narrow curvature range 
K < 1.2 for which p > 1, in accordance with their analysis assumption 
mentioned.- 'The non-zero slope of the Pg gl - K curve in the neighborhood 

* For, say, a/h = 600, K = 5 corresponds to panel central angle 

4 14°. 

o 
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of p = 1 and K « 1.2 suggests that P ssl < 1 for larger values 
of K. That this is the case is confirmed by the unconnected points 
shown in Fig. 3. These points are bifurcation bucling points 
that were taken from Pope's (Ref. 9) postbuckling curves for long 
and narr?^ SSI panels. Rehfield and Hallauer used these points to 
extrapolate their SSI results (see Ref. 8, Fig. 2). into the larger 
K region for which .p ss] _ < 1. Thus, Pope's limited number of SSI 

results, show, perhaps somewhat surprisingly, that the panel buckling 

■ l Ctfij. bljL ^ _ 

loadsV'^xcept^or^Th^reTativdly extremely 'narrow panels (K £ 1.2). 

The smallest value of the non-dimensional load parameter obtained 

from Pope's SSI results, is p^. .6 for K ^ 2.8.. However, the trend 

of Pope's results with increasing K suggests that an even lower 

value of p will be realized for a larger value of K. The existence 

of such low panel buckling loads is somewhat analogous to the well 

known low buckling loads for simply supported complete cylinders 

0 

with no circumferential constraint (v ^ ft, = 0, "weak in shear 

boundary condition) on the curved edges. For that case, the cylinder 
buckling load is approximately one-half (Refs. 10-12) of the 
classical cylinder buckling load (for which w “ M x = N x “ v 51 °> • 
However, for the complete cylinder, the "weak in shear" boundary 
condition rarely occur in practice (Ref. 13), whereas, in contrast, 
the SSI panel is of some practical interest since free "in-plane" 
movement of the straight edges can be simulated ^Experimentally . 

Also of interest for the SSI panels is whether the SSI panel load 
ever reaches the full cylinder buckling load for sufficiently large 
K with $ o ■< 2*. This will also be investigated herein. 
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1.3 SCOPE OR PRESENT INVESTIGATION ... 

In view of the above discussion, it seems desirable to per- 
form a comphrensive analysis that covers the complete range of 
panel widths (0 < <f> ' < 2tt) . This is the aim of the present paper. 
Buckling loads and mode shapes will be presented for panels with 
the eight sets of boundary conditions defined in Egs. (1) and (2). 

The panels will not necessarily be assumed to be long; buckling 
loads will be generated for different L/a values. The results will 
be based on both Donnell-type linear shell theory and a more complete 
shell theory. As was also employed in the previously cited re- 
ferences, a linear membrane analysis is used for the prebuckling 
solution and the effects of initial inf erfections are not considered. 
The "exact” results presented herein were obtained from the BUCLASP 
2 computer program (see Appendix A) . 
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2. NUMERICAL RESULTS AND PARAMETRIC STUDIES 


The present investigation is aimed at locating the criti- 
cal angle, l P cr r where <f> cr is the smallest value of $ for which 
the panel axial buckling load is identical with that of the 
complete cylinder. Studying of the influence of the different 
combinations of in-plane boundary conditions along the straight 
edges of the panel, namely SSI to SS4 and CCl to CC4 (Equations 
(1) & (2) } on the location of $> and the panel critical axial 

C3T- 

loads. Parametric studies on the effect of panel geometry, 

(L/a) and (a/h) on the critical load and verification of the 


2 2 1/2 2 

existance of a panel geometry parameter K = (1/tt ) [3 (1-v )] (b /ah) 


(1/tt ) [3 (1-v ) ) 


1/2 


(a/h) 4> , analog to the Batdorf shell curvature 


t 2 

parameter Z = (L /R^) , which has been defined in 


References 5 and 8. 

A main shell geometry, the so called "MARSHALL unstiffened 
Cylinder", has been chosen for the present thorough investigations. 
The dimensions and properties of this shell are as follows: 

L = 94”; a = 60"; h - ,1"; E = 10 7 p.«s.i. and = 1/3 (L/a - 1.567 
a/h - 600) 

For the parametric studies L/a has been changed to .1 and 5.0 and 
a/h to 100 and 2000, thus allowing studies with short and long 


panels as well as thick and thin ones. 



In Figures 4A and 4B, the ratio of panel buckling load over 
the complete cylinder buckling load, p = Nx o^^ x ^cyl (obtained ky 

BUCLASP for Donnell type analysis), versus the panel angle, $, 
is shown for the "MARSHALL" type panel. Figure 4 A shows the in- 
fluence of the SS type boundary conditions and Figure 4B that of 
the CC type boundary conditions. It is seen from these figures 
that the most effective in-plane restraint is X/ r =0 along the 
straight edges of the panel resulting in higher and almost identi- 
cal loads for the SS2 and SS4 as well as CC2 and CC4 boundary 
conditions. However, Figure 4B reveals that the CC curves are 
closer together than the SS curves in Figure 4 A which indicates 

4 * 

that the in-plane boundary conditions are more influential for 
the SS type boundary conditions. Also, the prevention of out 
of plane edge rotation dominates and hence is more important than 
the condition M'-O . 

It can also be observed in Figure 4B that for narrow panels, 
clamping of the straight edges results in higher critical loads 
than for simply supported edges. 

It is observed in Figures 4 b that in the case of CC boundary 
conditions, all of the curves approach p~l, the complete cylinder 
buckling load from above and the smallest value of 6 is observed 
for the CCl boundary conditions with <f> 10°. It is seen from 

this figure, that <i> varies in the range, 10 g <J> £.20 for all 

of the sets of in-plane boundary conditions. In contrast, 

Figure 4 a shows that in the case of SS boundary conditions only 
SS2 :to SS4 approach the complete cylinder critical load from 



above whereas the SSI boundary conditions are poorly behaved; 
cut the line p=l at about <j>*6.5°, decrease to a minimum of 
p^.465 for <{>325.0 and then increase to approach an asymptotic 
value of psr. 81'9 rather than p=l. This behavior is qualitatively 
similar to that experienced for "lightly" stringer-stiffened 
shells (Reference 14) with weak in-shear in-plane boundary con- 
ditions. Conclusively, the poorly behaved set SSI boundary 
conditions makes the CC1 set much more prefferable to the ex- 
perimentalist, but of course he would have to guarantee the pre- 
vailing of the CC1 boundary conditions. This also calls for the 
analysis of the combined effect of out of plane rotational springs 
together with the different sets of in-plane boundary conditions 
on the panel critical load. 

The critical loads were calculated with the aid of BUCLASP 
both for Donnell and Flugge type stability equations and the re- 
sults are presented in Table 1 of Appendix B* It can be ob- 
served in this Table that the results obtained by Flugge type 
equations are close to the ones discussed above, with the follow- 
ing exceptions : the CCl boundary conditions also reveal a ten- 

dency to approach the critical value from below, but with a value 
of p very close to p — 1 , only the SS3 boundary conditions ap- 
proach the value p=l from above whereas SS2 , SS4 and CC2 to CC4 
approach a value of p slightly above p— 1 and hence there is 
actually no existence of <f> cr for these sets of boundary conditions. 

<9 

It is also seen that for a panel with <{>=3 60, Donnell type aha- 
lysis predicts an axially multiwave buckling mode (excluding SSI) 
whereas Flugge type analysis predicts half a single wave mode. 



In Tabels 2 to 5 of Appendix B, the calculated results for 
a/h = 100; a/h = 2000, L/a = .1 and L/a = 5.0 are presented 
respectively. 

Plotting of the results for the "thick” panel (a/h) = 100 
in a similar manner to Figures 4 a and 4B show that the in-plane 
boundary conditions are less effective for this type of panels - 
almost no effect for the clamped boundary conditions and less 
influence of the SS3 boundary conditions when compared with the 
SS2 and SS4 boundary conditions. The SSI boundary conditions 
behave similary to the "MARSHALL” type panels. The discussion on 
the correlation of the Donnell and Flugge buckling loads applies 
also for this type of panels, except for the very narrow panels 
with SSI and CCl boundary conditions. It is found that the 
Flugge type analysis results in an in-plane Euler buckling load 
not included in Donnell type analysis. vHence the Flugge criti- 
cal loads are much lower. than the Donnell ones and it can be 
shown that when this mode becomes critical for a "thick" panel, 
the second Flugge buckling load corresponds to the Donnell 
critical load for the same panel. In the CCl case there are 
no "pure" Euler modes because clamping along the straight edges 
imposes the condition v=-w,^. 

The results for the thin panels, a/h - 2000 - Table '3 of 
Appendix B-reveal that these panels behave very similar to the 
"Marshall" type panels, abrupt curves obtained for this except 

for the SSI boundary conditions where an abrupt change in the 

» 

value of p is noticed for 4>=360. This value of p = .977 contra- 
dicts the value of p=.805 obtained by the Flugge type analysis. 



It is also found that for this type of panels, the SSI, SS3, 

CC1 and CC4 boundary conditions all approach p=l from below 
with a value of p very close to unity. 

Plotting the results for the short panels, L/a = .1 
(See Table 4 of Appendix B) results in conclusions similar to 
those obtained for the "MARSHALL" type panels, except for the 
SSI boundary conditions where p>l for <f> = 360° and p does not 
reach the minimum values obtained for the cases discussed above. 
Also, results obtained by Flugge type analysis are in better 
agreement with the Donnell type results for this panel configu- 
rations . 

From Table 5 of Appendix B it is found that., for long 
panels, L/a = 5.0 the behavior of the panels is also similar 
to that of the "MARSHALL" type panels. No correlation between 
Flugge and Donnell buckling loads has been obtained for the 5° 
panel and SSI boundary conditions. Calculations show that the 
Flugge critical load corresponds to an in-plane Euler buckling 
load, which Donnell's analysis excludes. It is also found that 
the CCl curve approaches p=l from below. 

It should also be noted in Tables 1 to 5 that the results 
obtained with aid of BUCLASP for SS3 and CCl boundary conditions 
are in excellent agreement with those obtained by the close form 
solution, Equations (4) and (9) . 

As stated previously the present paper is aimed at veri- 
fying the existance of the panel geometry parameter K of 
Equation (7). In Figures 5, p has been plotted versus K for 
SS boundary conditions and in Figures 6 for the CC boundary con- 



ditions . All of these figures except 5A r for the SSI boundary- 
conditions, indeed verify the existance of K independent of 
(a/h) and the panel angle, $ . For each set of boundary conditions 
a single curve is obtained. Note that the results for SS4 and 
CC4 boundary conditions are not included as they coincide with 
those for SS2 and SS4 boundary conditions respectively. 

A similar study has been performed on the length effect and 
is presented in Figures 7 and 8. These figures reveal that only 
for short panels, L/a = .1, such an effect exists. 

The representation of Figure 5A has not verified the exis- 
tance of the unique parameter K for the SSI boundary conditions. 
Hence, instead of presenting the results for this type of bound- 
ary conditions as p vs. K like in Figure 5A, an attempt has been 

. N, 

xo 

made to present the results in the form X = ^ — vs. K on a log 

P 1 

representation. , This is shown in Figure 9 and it is observed 
that by this kind of representation K does also become a single 
parameter for the SSI boundary conditions. 
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.APPENDIX A 


A- 1 BUCLASP 2 

The numerical results presented in this investigation were 
obtained from the BUCLASP 2 (Buckling of Laminated Stiffened Plates) 

A t ^ , , 

computer program . BUCLASP 2 is applicable to stiffened prismatic 
structures composed of composite flat plate, cylindrical panel, 
and beam " el e-men t s" . Classical simple support boundary conditions 
(Eq. (3)) are assumed for the curved edges, and a linear membrane 
prebuckling analysis is employed. As a consequence, the axial 
buckling modes are sinusoidal and a truly one-dimensional analysis 
is effected through the use of Fourier series representations for 
the axial variations of the buckling displacement components, u, v, 
w. Now since the structure is assumed to be prismatic, the equations 
governing the buckling behavior of each plate, panel, or beam com- 
ponent of the structure possess constant coefficients, and thus 
"exact" stiffness equations are readily obtained for each element* 
These elements are then assembled together through the use of the 
direct stiffness method to yield the "exact" stiffness equation for •• 
the entire structure. The elements referred to herein are the 
structural components that occur naturally in the stiffened structure, 
such as plates, panels, and beams, and they are not obtained through 
a spatial discretization as is done in the finite element approach. 
Buckling loads are then obtained upon specification of the boundary 
conditions along the straight edges of the structure. 


"Elastic Buckling Analysis for Composite Stiffened Panels and 
other Structures Subjected to Biaxial Inplane Loads" , by A. V. 
Vis war. a than, M. Tamekuni, NASA CR-2216, March 197 3. 
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A. 2 GOVERNING EQUATIONS 

Specialization of the more general equations presented in 
BUCLASP 2 to the case of the unstiffened, isotropic, cylindrical 
panel structure considered herein results in the following con- 
ventional Donnell shell equations (see Figs. 1 and 2 for sign 
conventions): . 

EQUILIBRIUM 


N + N = 0 

x,x xy,y 


N + N = 0 
y,y 


. N 
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+ 2M + M 


y _ 


N W 


x,xxx xy,xy y ,yy a xo ,xx 


CONSTITUITIVE 


N x = (£ x + ve y ) 

1-v 


M x = E ^ ' T~ (K x + ' ,K y ) 

x 12 (1-v ) X y 
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y y x 
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y 12 (i-v 2 ) y x 


N “ Ghy 
xy 1 xy 


M 


Gh' 


xy 


12 (1-v ) 


: — n. 

, xy 


12 



KINEMATICS- 


e x = U 'x 


= - W 


W 

e = v - — 
y rY a 


K = - W 

y . »yy 


Y = U + V K — — 

xy ,y ,X xy 


Results are also presented herein based on a sow «s?lete shell 
theory that is presented in the 3 U CLASP 2 
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ABSTRACT 


Boron infiltration of extruded aluminum stiffeners is a selective 
reinforcement concept being studied at the Langley Research Center for 
specialized aerospace applications including an intertank skirt and a wing 
box panel. One critical design problem associated with this concept is the 
transfer of loads from end attachments, joints, or splices into the boron 
reinforced stiffener. The abrupt increase in stiffness in end regions causes 
high shear stresses to develop: (l) in the metallic weh member connecting 

the end attachment and the boron reinforcement, and (2) in the adhesive bond 
between the infiltrated boron and aluminum stiffener. High shear stresses in 
the end region, for example, contributed to premature failure of one previ- 
ously tested boron infiltrated panel. 

To reduce critically high web shear stresses, tapering is proposed so 
that the axial stiffness of the boron reinforced rod is gradually increased 
from zero stiffness at the free end. This is physically accomplished in a 
post- infiltration machining or grinding operation. Critically high epoxy 
bond shear stresses can be reduced by a grooving operation which removes boron 
material but retains most of the bond surface area. A linear taper used to 
reduce web shear stresses, unfortunately, increases the amplitude of bond 
shear stresses, and a combination of tapering plus grooving may be n.ecessary 
to satisfy both web and bond allowable shear stress constraints, 

This paper presents results of analytical and experimental studies con- 
ducted to evaluate the effectiveness of end tapering, grooving, and end 
fixture sculpturing to reduce critical shear stresses in boron infiltrated 
extruded structures. The two configurations studied were a bi- element tension 
specimen and a wing box compression panel. 



LOAD INTRODUCTION TECHNIQUES FOR BORON INFILTRATED ALUMINUM PANELS 
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and 
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Cincinnati, Ohio 


INTRODUCTION 

NASA has studied, through contract (Ref. 1) and in-house efforts, the 
selective reinforcement concept in which collimated boron and epoxy resin are 
infiltrated into cylindrical voids in extruded aluminum alloy stiffeners (see 
typical stiffener cross sections in Fig. i). The primary advantage of infil- 
tration compared with bonding high modulus plies directly onto the surface of 
a metallic stiffener is the protection provided the boron/epoxy by the alumi- 
num annulus. To evaluate the infiltrated stiffener concept numerous strength 
and crippling specimens have been tested (Ref. 1) and several large compo- 
nents are in various stages of being built and tested. Aluminum stiffeners 
have been successfully extruded in 6005 and 7075 series alloys. Cross- 
sectional conf igutations have included hats, channels, and "T" and "Y" 
sections. 

It is necessary for most practical structural applications of the boron/ 
epoxy infiltrated stiffener concept to incorporate a metallic load introduc- 
tion end fixture. Loads are carried totally by metal at the end of the, 
stiffener and are transferred into the boron through shear in the metallic 
veb of the fixture and stiffener. The maximum length which boron /epoxy may 
be infiltrated into the annulus of a stiffener has not been established; how- 
ever, some applications may employ splice joints which also require loads to 
be transferred into the boron from an all-metal section. Since the boron 
typically carries 50 percent or more of the total load, attention must be 
given the design of the load introduction fixture to prevent premature fail- 
ure and to minimize the fixture weight. 

The present paper summarizes a study of the load transfer problem. The 
study was prompted by evidence in early test specimens of premature failure 
in the load introduction region. Recent analytical studies reveal unaccept- 
ably high shear strsses in the load introduction region of a proposed wing 
box panel application. One way found to reduce high shear stresses is to 
reduce the axial stiffness of the boron/epoxy infiltrated aluminum annulus 
near the stiffener free end. Boron/epoxy and aluminum material are removed 
from the stiffener either in a tapered grinding or a grooving operation. 

This paper reviews the performance of infiltrated components, presents 
the philosophy for tapering and grooving to reduce end region shear stresses, 
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describes an analytical and experimental investigation of a double annulus 
tension specimen to assess the merits of tapering, and finally describes how 
tapering reduces the shear stresses to acceptable magnitudes in the wing box 
panel design. 


SYMBOLS 

Length 

Axial load per unit width 
Axial dimension 
Axial stress in boron 
Maximum axial stress in boron 


Aluminum web shear stress 
Maximum aluminum web shear stress 


Epoxy bond shear stress 
Maximum epoxy bond shear stress 


PERFORMANCE OF INFILTRATED COMPONENTS 

Photographs of several boron infiltrated stiffened structural components 
are presented in Figure 1. The cylindrical interstage tank specimen shown in 
Figure 1(a) is I .96 m (77 in.) long and has a diameter of 3.91 m (154 in.). 
External longitudinal stiffeners are hat sections with boron infiltrated into 
all four corners. Internal ring stiffeners have lr H n shaped cross sections 
with boron located only in the top half of the stiffener. See inserts in 
Figure 1(a). The structure is moderately loaded with an axial compressive 
load equal to 823 kN/m (4700 lbf/in. ). This 'shell structural specimen has 
been recently tested successfully in the 3*39-MmN (30-million in. -lb) bending 
test fixture in the Structures Laboratory at Langley Research Center. The 
DC-10 floor strut (Fig. 9(b)) is being evaluated under a joint AVCO 
Corporation/McDonnell- Douglas Corporation program and is reported upon in 
another paper in these proceedings. 

The two compression panels shown in Figures l(c) and 1(d) are constructed 
with boron infiltrated m Y" stiffeners. These panels weigh 1.18 kg/m 2 
(2,6 lbm/ft^) and are designed to carry an axial load of 1.26l MN/m 
(7200 lbf/in. ). The panels are 1.22 m (48 in.) long and 0.86 m (34 in.) wide 
and differ only in the end fixture design. The compression panel shown in 
Figure 1(c) has massive steel end fixtures bolted to the skin of the panel 
and to the web of the "y" stiffners. Loads are transferred into the 
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boron/epoxy through shear flow In the metallic stiffener web. n 

end support simulates a simple support end boundary condition. The specimen 
was tested at Langley Research Center and failed at approximately 65 percent 
of the design ultimate lead. Premature failure was caused in part by high 
.w stresses in the load introduction region of the stiffeners. 

t closeup view of the shear failure in the end region is shown in the 

photograph in Figure 2. Two characteristic types of shear Jfen^^as^ailed 
seen in the photograph. First, the aluminum web of one stiffener 
in shear as evidenced by the physical separation of the top of the stif 
from the web (a coin has been slipped between the two separ ated pa ^ Se=" 
ond, the epoxy bond between the boron and the aluminum has failed on another 
stiffener allowing the boron rod to displace relative to the end of th 
aluminum extrusion. . 

The end fixture for the second compression panel is designed to permit 
loads to be transferred directly into the boron/epoxy (Fig. 1 (d)). This tlat- 
end test condition permits assessment of the load-carrying ^pabilityofthe 
panel without the complications imposed by the load transfer prob p 

stiffeners. This specimen was also tested at Langley and failed in an Euler 
mode at 91 percent of the design ultimate load. No shear type failures in 
stiffeners were observed. 

The latest effort in boron infiltrated structures is a wing box cover 
panel which is currently in the late stages of desipi. The s ,^ c ® /. y 

designed to be heavily loaded (axial loading equal 4-90 Mh/m (..8,000 M/m. ) 
and has a sculptured load introduction fixture represer Native °t th ; ^ ^ e |uiie. 
for a ioint in the wing box application. The assembly process for this struc 
t^e involves T^gsten^nert Gas (TIG) welding a T" stiffener to the web of 
an integrally machined aluminum plank. Shear stresses m the web weld regio 
are especially important because the weld process reduces the allowable . 
strength of the 2219 aluminum alloy. Initial analytical studies of the wing 
box ifad introduction region indicated shear stress magnitudes in the web 
did exceed the ultimate shear strength. The end tapering approach for 
reducing web shear stresses was applied to this configuration and resu ,s aie 

reported in subsequent sections. 


END TAPER CONCEPT DESCRIPTION 


The shear stress problem associated with the transfer of loads into the 
concentrated mass of boron/epoxy in an infiltrated stiffener is basically a 
twofold problem. First, the abrupt increase in stiffness in the end region 
causes high shear stresses to develop in the metallic web member connecting 
the end fixture and the infiltrated boron reinforcement. Second, high shear- 
stresses can develop in the epoxy bond interface between the aluminum, and the 
boron /epoxy materials. 


Drawings in Figure 3 show schematically (with distortions greatly exag- 
gerated) how boron stiffener end tapering reduces the shear stress magnitude 
in the connecting metallic web member. When the encapsulated boron is 



terminated by a square cut, loads carried by the skin are transferred into 
the boron over a .very short length near the end region, thereby causing high 
shear stresses to develop in the aluminum web. Straight lines drawn on the 
web of the deformed structure represent the displaced position of lines which 
on the undeformed structure would originally be parallel and uniformly spaced. 
Large rotations of these lines between undeformed and loaded conditions 
represent high shear stress regions. The shear stresses are largest near the 
free end and decrease to zero at a point along the stiffener axis where the 
load redistribution is complete. 

In the tapered approach, illustrated at the bottom of Figure 3* the 
stiffness of the encapsulated boron/epoxy rod is gradually increased from 
zero at the free end to the nominal value at the termination of the taper. 

This arrangement result's in nearly uniform shear stresses in the web end 
region (represented by rotated nearly parallel lines in the deformed struc- 
ture), and the magnitude of stress for the tapered structure is less th. .n the 
maximum shear stress in the web of the square- cut terminated stiffener. This 
gradual buildup of stiffness is physically implemented by machining or grind- 
ing a taper in the boron infiltrated rod to remove material for a few centi- 
meters near the end. In the present effort a linear taper is employed. A 
typical example of a linear taper is shown on the top of the specimen shown 
in the photograph in Figure 4. The tapering approach presented here for the 
infiltrated boron concept is basically analogous to the approach proposed by 
deBruyne (Ref. 2 ) in 1944 for bonded or glued lap joints. deBruyne found that 
tapering the two ends of a lap joint doubles the failure load compared to 
nontaper ed joints. 

As previously mentioned, the shear stress in the epoxy bond interface 
may also be critical. Critically high bond shear stresses can be reduced by 
a grooving operation such as that shown on the bottom of the specimen in the 
photograph in Figure 4. Grooving reduces the rate at which load is trans- 
ferred into the boron while retaining most of the bond surface area. 

A linear taper which is effective in reducing web shear stresses also 
(unfortunately) increases the bond shear stress. This increase occurs since 
stress is defined as force per unit area and the bond area lost through taper- 
ing is reduced by a greater factor than the corresponding rate at, which load 
is transferred through the bond into the boron. For structural applications 
in which both web and bond shear stresses are critically high, one possible 
solution is to machine a linear taper followed by a grooving operation which 
removes boron but leaves the bond surface area intact. 


TENSION SPECIMEN 

Specimen and Experiment Description 

A tension specimen used to study the load transfer problem is shown in 
the photographs of Figure 5* The specimen consisted of a 6005-T5 aluminum 



bi- element extrusion which had an overall length of 48.3 cm (19 in. )• The 
web was removed for 10.2 cm (4 in.) in the center of the specimen resulting 
in a load introduction length of 19.1 cm (7*5 in.) at either end. This cut- 
out permits loads introduced at the end tabs to be fully transferred into the 
boron reinforced stiffener in the specimen center. Each aluminum annulus was 
filled with 857* 0.14 mm { 5 . 6 mil) diameter boron filaments and subsequently 
infiltrated with a room- temperature curing epoxy resin. Steel load introduc- 
tion straps 2.34 mm (0.1 in.) thick were bonded and bolted to the web of the 
test specimen in a symmetrical double lap configuration. The cross section 
of the two boron infiltrated stiffeners was reduced at one end of the specimen 
by grinding a linear taper 9*54 cm (3-75 in.) long. The boron at the other 
end of the specimen was terminated by a transverse square cut. This configu- 
ration permitted comparison on the same specimen of the load transfer response 
for a tapered and nontaper ed end. 

Geometric constraints make it impractical to measure experimentally the 
shear strains in the aluminum web and epoxy bond. Instead; axial strains 
were measured on the surface of the boron stiffened aluminum annulus using 
38 strain gages positioned along the specimen length. Comparison was made 
between experimental and analytical results and the analysis was then used 
to determine shear stresses in the aluminum web and epoxy bond. Bond shear 
stresses were determined by calculating the rate at which loads are trans- 
ferred into the boron and dividing by the circumferential bond dimension. It 
was assumed that the bond shear stress is uniform around the circumference. 
Tests were performed using a 445-kN (100,000- lb) capacity hydraulic testing 
machine and strains were recorded using an automatic data acquisition system. 


NASTRAN Model 

The finite- element computer program NASTRAN (NAsa STRuctural Analysis) 
was used to analyze the shear stress load introduction problem. Loads which 
stress the boron. to its ultimate strength may result in stress concentrations 
which produce plasticity effects in the aluminum. In the current study, only 
a linear analysis was conducted; however, plasticity effects usually reduce 
peak stresses so that an elastic analysis is conservative. The study was 
restricted to mechanical loads; stresses due to potential thermal expansion 
differences between the aluminum and boron (due to curing or temperature 
effects) were not considered. 

A schematic of the finite- element model used to represent the boron 
reinforced bi- element specimen is shown in Figure 6 . Rectangular constant 
strain elastic membrane plate elements were used to represent the aluminum, 
boron, and steel components. Although the full model is shown in Figure 6 , 
symmetry about three axes was used to reduce the computational model size. 

A total of 264 plate elements and 191 grid points were used in the computa- 
tional model. No attempt was made to represent details of the aluminum 
annular cross-section geometry. The aluminum and boron materials were assumed 
to have moduli of elasticity of 68 . 9 GPa (10 million psi) and 206.8 GPa 
(30 million psi), respectively. A uniform tension load was imposed on the 
steel straps at the specimen ends. 



Analytical and Experimental Results 


Comparison of analytical and experimental results for the axial stress 
in the boron rod of the bi-element tension specimen is presented in Figure 7« 
Results are presented for both tapered and nontapered cases. Experimental 
results shown were calculated using strain-gage data by assuming the boron and 
aluminum strains are identical. The stress magnitude has been normalized by 
the maximum stress amplitude (a - 800 MPa ( 116,000 psi)) which occurs at 

max 

the center of the specimen where the load is completely carried by the infil- 
trated boron circular stiffener. Results are presented for the complete 
specimen to establish correlation between experiment and theory. Shear stress 
studies, however, are focused on the end region (i.e v , x/L < 0.5) since this 
region is representative of the load introduction problem. 

Correlation between experimental and analytical results is reasonably 

good. The effect of tapering in the end region is to increase the stress in 

the boron over the nontapered result (Fig. 7)* 'This increase, as explained 

earlier, occurs since the bond area lost through tapering is reduced by a 

greater factor than the corresponding rate at which load is transferred 

through the bond into the boron. The total force in the boron and the rate 

at which load is transferred into the boron, however, is reduced for the 

tapered case. This situation results in a reduction in the maximum, shear 

stress in the aluminum web for the tapered end of approximately 60 percent ' 

as can be seen in the NASTRAN results presented in Figure 8. When the boron 

rod at the center of the test specimen is stressed to = 1.58 GP& 

(200,000 psi), the maximum web shear stress for the nontapered case is 

o = 128 MPa (18,500 psi). The ultimate shear stress of 6005 aluminum is 

xy 
max 

158 MPa (20,000 psi) (Ref. 5) and the web shear stress, therefore, is not 
critical at this load even for the nontapered case. 


The preceding example demonstrates the effectiveness of boron stiffener 
tapering to reduce high web shear stress in the end region of boron infil- 
trated structures. A comparison of the epoxy bond shear stress in the end 
region for the tapered and nontapered tension specimen is presented in Fig- 
ure 9* As indicated previously, the net effect of tapering is bo increase 
the maximum bond shear stress. In this example with the boron stressed at 


0 s= 1,38 GPa (200,000 psi) in the center test section, the tapered end maxi- 
A = 19,5 MPa (2800 psi). This value is 


mum bond shear stress is t 


xy. 


max 


approximately 100 percent higher than the maximum bond shear stress for the 
nontapered end. The recommended allowable shear stress for the epoxy bond' 
used in these specimens is approximately. 15- 2 MPa (2200 psi). Tests conducted 
by the AVC0 Corporation, however, indicate shear stress concentrations as high 
as i-l.lj- MPa (60O0.psi) can be carried by the bond. 


The analytical assumption of constant shear stress in the bond is non- 
conservative. The maximum bond shear stress near the intersection of the 
stiffener with the web is greater than the average bond shear stress and is 
reduced in magnitude from that point around the bond circumference. If a 



more refined analysis were conducted which included local bond shear stress 
gradients, the maximum bond shear stress would increase by a greater amount 
for the nontapered than the tapered case and, in effect, would decrease the 
relative disadvantage shown for a linear taper. 


WING BOX PANEL 
Structure Description 

The wing box panel structure consists of boron reinforced extruded 
stiffeners which are TIG welded to an integrally machined aluminum plank. 

The end fixture is also machined integrally with the wing plank and is 
sculptured to permit the gradual transfer of loads into the boron reinforced 
stiffener. The concept and weld location is shown schematically in Figure 10. 
The weld is made by joining two extrusion elements to the web of the inte- 
grally machined plank to form a M T" stiffener. The compression panel of the 
wing box is designed to carry an ultimate axial compression load of 4.90 MN/m 
(28,000 lbf/in. ). Constraints imposed by welding, heavy loading, and a 
requirement to show a weight savings over an all-metal design combine to make 
the load introduction an important design problem. The weld area is critical 
since TIG welding reduces the allowable shear stress for the 2219 aluminum 
material used in the wing plank from 262 MPa (38, 000 psi) to about 138 MPa 
(20,000 psi). 

Analyses of the wing box panel preliminary design showed the shear stress 
in the vertical web weld region to be substantially greater than 13 8 MPa 
(20,000 psi) for a boron reinforced stiffener terminated by. a square cut. End 
region tapering was studied to determine the capability of this approach to 
solve the problem. 


NASTRAN Model 

Drawings of the NASTRAN model used to represent the load introduction 
region of a typical wing box stiffener are presented in Figure 11. Selected 
cross sections show the sculpturing of the load introduction fixture. Sym- 
metry about the mid.plane of the "T“ stiffener was utilized to reduce the 
computational model size. Three-dimensional constant strain solid elements 
were used to model the thick end fixture and *skin components. A total of 
727 constant strain plate and solid elements were used in the computational 
model to represent the vertical and horizontal webs and 36 bar elements having 
axial stiffness only were used to represent the boron infiltrated circular 
stiffener. 

The wing box. compression test panel is '2.44 m (96 in. ) long, 0.91 m 
(36 in.) wide and has rib supports every 0. 6l m (24 in.). Although the load 
introduction model studied was only 0.49 m (19-4 in.) long, analyses indicate 
that most of the load transfer takes place within this length. In the analy- 
sis, loads were assumed to be applied by imposing a uniform displacement to 
the end of the load introduction fixture and restraining the axial 



displacements at the other end of the model* The aluminum and boron were 
assumed to have moduli of elasticity of 68,9 GPa (10 million psi) and 241 GPa 
(55 million psi), respectively. 


NASTRAN Results 

The effectiveness of end tapering to reduce the shear stress in the 
vertical web weld region for the wing box compression panel loaded with an 
axial load of 4.90 MN/m (28,000 Ib/in. ) is presented in Figure 12. The maxi- 
mum shear stress a for the nontapered case of 207 MPa (30,000 psi) was 

^max 

used to normalize the ordinate and the analytical model length L of *4-9*3 cm 
(19.4 in.) was used to normalize the abscissa. The shear stress for the non- 
tapered case exceeds the 13 8 MPa (20,000 psi) allowable shear stress by 
50 percent. Linear tapers of 0.171 1 (( 8 , *4-1 cm (3*31 in.)) and 0.325 L 
(l 6.0 cm ( 6.31 in.)) reduce the maximum shear stress in the weld region to 
0.9 a (186 MPa (27,000 psi)) and O .67 a (138 MPa (20,000 psi)), 

^max 

respectively. The 0-325 L taper permits loads to be transferred into the 
boron without exceeding the linear elastic shear stress allowable. The peak 
shear stress near the end for the nontapered case has been replaced for the 
0.325 L taper by a nearly uniform stress for the first 0.21 L ( 10 , 2 cm 
(4 in.)) of the stiffener. 

The epoxy bond shear stress for the 0.32-5 L taper and nontapered cases 
for the wing box panel loaded with an axial compressive load of 4.90 MN/m 
(28,000 lb/in. ) is presented in Figure 13* The maximum shear stress for the 
0.325 L linearly tapered case is 29*6 MPa (4300 psi) and this value has been 
used as t . The recommended allowable bond shear stress value of 
■^max 

15.2 MPa (2200 psi) (0.51 r ) is not exceeded for the nontapered case. 

^max 

The 0.325 L taper increased the bond shear stress compared to the nontapered 
result by approximately 53 percent and the value exceeds the recommended 
allowable by approximately 50 percent. 

The effectivennss of a combination taper and groove to reduce the bond 
shear stress was studied analytically by adding a groove to the 0.325 h 
taper case. The groove was assumed to be rectangular in cross section and 
3.8 mm (0.15 in.) wide. The groove was begun at a distance 0. 077 L (3«8l cm 
(l .5 in.)) from the stiffener end and continued parallel to the stiffener 
over a distance of 0.103 L (5*1 cm (2 in.)). The groove slope from that 
point was such that the bottom of the groove intersected the termination point 
of the taper. The reduction in bond shear stress which results is presented 
in Figure 13 . A substantial reduction occurs in the region of the taper. 

The groove cross section was insufficient, however, to reduce the bond shear 
stress to the 15.2 MPa (2200 psi) recommended allowable. To be fully effec- 
tive, the taper should be initiated around 0.036 L (1.8 cm ( 0.7 in*)) from 
the stiffener end near the point where the tapered boron begins and additional 
boron should be removed. 



The small change in stiffness provided hy the groove modifies only 
slightly the shear stress distribution in the aluminum web weld region. This 
result can be seen in the comparison of tapered and tapered plus grooved web 
shear stress distributions plotted in Figure 12. 


DESIGN IMPLICATIONS 

Most practical structural applications of boron Infiltrated stiffener 
selective reinforcement concepts require an end fitting in which loads are 
transferred through shear from an all-metal Joint into a cross section con- 
taining boron reinforcement. High web shear stresses can be reduced by taper- 
ing the boron infiltrated stiffener cross section and high bond shear stresses 
can be reduced by grooving operation which removes boron but retains most of 
the bond surface area. A load introduction fixture design which meets both 
web and bond allowable shear stress requirements may require both tapering 
and grooving. Each structural application has unique design constraints and 
determination of a satisfactory load introduction fixture design may require 
several iterations. Note, however, that load introduction shear stress prob- 
lems can be solved by removing rather than adding material. 

The web which connects the skin to tfie boron reinforced annulus must be 
sufficiently thick (1) to carry the shear loads imposed in the load intro- 
duction region, and (2) to carry axial loads without crippling locally. For 
the wing box panel, the thickness required to meet shear stress requirements 
exceeded the crippling requirement. Constraints imposed by TIG welding 
restricted thickening the web in the end region and additionally reduced the 
allowable shear stress in this critical section. Alternate assembly tech- 
niques such as riveting also present design problems, especially for heavily 
loaded structural applications. This illustration emphasizes that the com- 
bined effect must be considered when making design decisions. 

The feasibility of both the taper and groove machining operations has 
been successfully demonstrated in the machine shop. Grinding a linear taper 
is a simple operation. Cutting a groove in a tapered stiffener requires 
greater precision since it is necessary to center the cut in order to retain 
the desired quantity of boron and leave the bond, interface surface undamaged. 
Accessibility to the boron stiffener free end must be provided when "tight 
fit" geometric constraints are imposed such as in the wing box panel load 
introduction fixture. 

Shear stress problems can also develop in the end region from -thermal 
expansion incompatibility between the boron and aluminum. Typically, the 
contributions of mechanical and thermal loads superpose to amplify the prob- 
lem. The removal of web material in the vicinity of the free end is suggested 
in Reference 4 as a technique for separating the maximum bond shear stress 
amplitudes for these two types of loading. Web removal transfers the maximum 
bond shear stresses caused by mechanical loading away from the free end. 

Based on bond shear stress calculations (Ref. 4), the maximum diameter 
which can be used for a boron infiltrated stiffener subjected to a 222 K 
(4-00° F) temperature differential has been determined to be 0.71 cm (0.28 in.). 



For a room- temperature cured epoxy such as used in the current investigation, 
typical thermal excursions would be less than 222 K (400° F), Only limited 
work has been done on thermal cycling for the boron infiltrated stiffener 
concept, and more study is needed. 


CONCISIONS 

An analytical and experimental study has been made of the load transfer 
mechanics cf structures stiffened by boron infiltrated extrusions where loads 
must be transferred through shear in an aluminum web from a load introduction 
fixture into a boron infiltrated stiffener. It was found that critically 
high shear stresses in the aluminum vjeb can develop in the load introduction 
region, especially for heavily loaded structures. Removing material to 
reduce the axial stiffness of the boron infiltrated aluminum annulus for 
several inches near the stiffener end using a linear taper was found effective 
in reducing the aluminum web shear stress. Reductions in the web maximum 
shear stress of 60 and 30 percent, respectively, were demonstrated for a ten- 
sion test specimen and projected for a wing box panel. Analytical and experi- 
mental results for the axial stress in the boron for the tension specimen 
showed reasonably good agreement. A 16. 0-cm (6.31-in.) long linear taper was 
sufficient to allow critically high shear stresses in the aluminum web weld 
region of the wing box panel to be reduced to an acceptable value. 

Critically high shear stress in the bond between the boron /epoxy and the 
aluminum can be reduced by grooving a cut in the circular stiffener. Grooving 
reduces the shear stress in the bond by reducing the magnitude of force which 
must be transferred across the bond while maintaining a high percent of 
effective bond area. Although a linear taper reduces shear stress concentra- 
tions in the aluminum web, it increases the shear stress magnitude in the 
epoxy bond joining the boron to the aluminum annulus. The simultaneous 
reduction of web and bond shear stress magnitudes can be accomplished by a 
combination of tapering and grooving. 
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Pig-are 1. Boron infiltrated stiffener structural applications. 





(a) Nontape red. 



Figure 3. Tapering of "boron infiltrated stiffener to reduce critical web shear stresses. 
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Figure 5» Boron infiltrated bx-element tension specimen. 
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Figure 6. HASTRAN model of di-element tension specimen. 







Figure 7. Comparison of axial, stress in 
and linearly tapered end regions. L = 


boron of bi-element tension specimen for nontane^ed 
24.13 cm (9.5 in. ). a x = 800 MPa (116,000 pal). 

max 




Figure 8. Comparison of calculated shear stress in aluminum web adjacent to stiffener of 
bi-element tension specimen for tapered and nontapered end regions. L = 24.13 cm (9.5 in.) 

- 128 (18,500 psi) when - 0 in boron at specimen center equal I.38 GFa' 

^max x 

(200,000 psi). 



L 

Figure 9. Comparison of calculated shear stress in epoxy bond of bi-element tension specimen 
for tapered and nontaper ed end regions. L = 2k. 13 cm (9*5 in.), t = 19-3 MPa (2800 psi) 

'■’max 

when o x in boron at specimen center equal 1.J8 GPa (200,000 psi). 



(a) W: 

Figure 10. Wing "box integrally i 



(b) Typical stiffener. 
Figure 10. Concluded. 



(a) Assembled half model. 

Figure 11. NASTEAN model of sculptured load introduction region of wing 
box compression panel typical stiffener. 
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(b) Cross-section detail. 


Figure 11. Concluded. 
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Figure 12. Comparison of calculated shear stress in wing box weh weld region for tapered and 
nontapered stiffener end regions. L = 49*5 cm (19-4 in.), a — = 207 MPa (30,000 psi) 
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when specimen loaded with axial load of 4.90 MN/m (28,000 Ibf/in. ). 
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loaded with axial load of *4-. 90 MN/^ ( 28,000 Ibf /in. ) . 


